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Abstract 

In this paper, wc investigate a numerical algorithm for the pricing of swing options, relying 
on the so-called optimal quantization method. 

The numerical procedure is described in details and numerous simulations are provided to 
assert its efficiency. In particular, we carry out a comparison with the Longstaff-Schwartz 
algorithm. 
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qh; Introduction 



In increasingly deregulated energy markets, swing options arise as powerful tools for modeling 
^ ■ supply contracts [H]. In such an agreement between a buyer and a seller, the buyer always has 
. to pay some amount even if the service or product is not delivered. Therefore, the buyer has to 
manage his contract by constantly swinging for one state to the other, requiring delivery or not. 
P\| ■ This is the kind of agreement that usually links an energy producer to a trader. Numerous other 
1/-^ . examples of energy contracts can be modeled as swing options. Prom storages [51 [H] to electricity 
O I supply [T71 [7], this kind of financial device is now widely used. And it has to be noticed that its 
field of application has recently been extended to the IT domain |12j . 

Nevertheless, the pricing of swings remains a real challenge. Closely related to a multiple 
stopping problem |101 [9]. swing options require the use of high level numerical schemes. Moreover, 
^ ' the high dimensionality of the underlying price processes and the various constraints to be integrated 
■ in the model of contracts based on physical assets such as storages or gas fired power plants increase 
the difficulty of the problem. 

Thus, the most recent technics of mathematical finance have been applied in this context; from 
trees to Least Squares Monte Carlo based methodology [25l [lU [l8] , finite elements [26] and duality 
approximation [20] . But none of these algorithms gives a totally satisfying solution to the valuation 
and sensitivity analysis of swing contracts. 

The aim of this paper is then to introduce and study a recent pricing method that seems very 
well suited to the question. Optimal Vector Quantization has yet been successfully applied to the 
valuation of multi-asset American Options [21 [H [3] . It turns out that this numerical technique is 
also very efficient in taking into account the physical constraints of swing contracts. For sake of 
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simplicity we shall focus on gas supply contracts. After a brief presentation of such agreements 
and some background on Optimal Quantization methods [22], we show that a careful examination 
of the properties of the underlying price process can dramatically improve the efficiency of the 
procedure, as illustrated by several numerical examples. 

The paper is organized as follows: in the first section, we describe in details the technical 
features of the supply contracts (with firm or penalized constraints) with an emphasis on the 
features of interest in view of a numerical implementation: canonical decomposition and normal 
form, backward dynamic programming of the resulting stochastic control problem, existence of 
bang-bang strategies for some appropriate sets of local and global purchased volume constraints. 
Section 2 is devoted to some background on optimal vector quantization. In Section 3, our algorithm 
is briefly analyzed and the a priori error bound established in the companion paper [4] is stated (as 
well as the resulting convergence result of the quantized premium toward the true one) . In Section 4, 
numerous simulations are carried out and the quantization method is extensively compared to the 
well-known least squares regression algorithm "a la Longstaff-Schwartz" . An annex explains in 
details how the price processes we consider in this paper can be quantized in the most efficient way. 

1 Introduction to swing options 
1.1 Description of the contract 

A typical example of swing option is an energy (usually gas or electricity) supply contract with 
optional clauses on price and volume. The owner of such a contract is allowed to purchase some 
amount of energy qt^. at time tk, k = 0, . . . ,n — 1 until the contract maturity t„ = T, usually one 
year. The purchase price called strike price may be constant or indexed to past values of crude 
oil. Throughout the paper we will consider that the strike prices are constant and equal to K 
over the term of the contract. The volume of gas qt,. purchased at time tj is subject to the local 
constraint 

Qmin ^ 'Ztfc ^ Qrnax • 

The cumulative volume purchased prior to time tk{i.e. up to tk-i) is defined by Qt^. = YI^Zq Qtf It 
must satisfy the following global constraint (at maturity): 

n-l 

Qt = Qtk £ [Qmin, Qmax]- 
k=0 

Two approaches can be considered: 

- The constraints on the global purchased volumes are firm. 

- A penalty is applied if the constraints are not satisfied. 

The price at time t of the forward contract delivered at time T is denoted by Ft^r, {Fo,tk)o<k<n 
being a deterministic process (the future prices at time 0) available and tradable on the market. 

Let {St^)o<k<n be the underlying Markov price process defined on a probability space {^},A,F). 
Note that it can be the observation at time t^,, A; = 0, . . . , n of a continuous time process. Ideally 
St should be the spot price process of the gas i.e. St = Ft^t- However it does not correspond to a 
tradable instrument which leads to consider in practice the day-ahead contract Ft^t+i- 

We consider its (augmented) natural filtration J^^ = {^t^)o<k<n- The decision sequence 
(Q'tfe)o<fc<n-i is defined on {Q,A,F) as well and is .7^'^-adapted, i.e qt,, is Tt^, = (^{^tQ, . . . , St^.) 
measurable, k = 0, . . . ,n. At time the owner of the contract gets qt,,{St^. — K). 
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Remark 1.1. The results of this paper can also be applied to every physical asset or contract where 
the owner reward for a decision qt,, is a function tp{tk, qt^, St,^). In the case of supply contracts, 
i/j{tk, qtk, St^) = qt^{Stf. — K). As for a storage, qt^ represents the amount of gas the owner of the 
contract decides to inject or withdraw and the profit at each date is then 



k J 



-QtkiStk+Ci) if qtk>0 (Injection) 

-qt,^{St^ - cw) if Qtk<0 (Withdrawal) 

^/ Qtk = (Same level in the storage) 



where c/ (resp. cw) denotes the injection (resp. withdrawal) cost [6]. 

1.1.1 Case with penalties 

We first consider that the penalties are applied at time T if the terminal constraint is violated. For 
a given consumption strategy (qtk)o<k<n, the price is given by at time 

P(0, 5o, 0) = E ( 1^ e-^'^qt.iSt, - K) + e-''^P^(ST, Qr)|^o ) 

\fc=0 / 

where r is the interest rate. The function (x, Q) ^ Pt(^i Q) is the penalization: Pj,(x, Q) < and 
Pj,(x, Q)\ represents the sum that the buyer has to pay if global purchased volume constraints, say 
Qmin and Qmax, are violated. [6] have already investigated this kind of contract. 

Then for every non negative ^t^-i measurable random variable Qt^^ (representing the cumulated 
purchased volume up to tfc-i), the price of the contract at time tk,k = 0, . . . ,n — 1, is given by 

P(tk,St„Qt,) = esssup E ( |^ e-^(*^-**)%,(5i, - i^) + e-^^^-^'-^^P^C^T, gr)|5t J ■ (1) 
The standard penalization function is as follows: 

P^(X, Q) = -(Ax(Q - Qmin)~ +Bx(Q- Q„,ax) + ) (2) 

where A and B are large enough - often equal - positive real constants. 

1.1.2 Case with firm constraints 

If we consider that constraints cannot be violated, then for every non negative ^tk-i measurable 
random variable Qt,. defined on (Tl, A, P), the price of the contract at time i^j ^ = 0, . . . , n — 1 is 
given by: 

P(tk,St„Qt,)= esssup ^ E(^|^e-^(*^-*^)gt,(5t, -K)|5i,y (3) 

where 

j^Qmzn,Qmax _ I (gij J , ; (Q^J^^^^F) [qmimQmax], 9t£ ^ [(Qmin — Q)^ ,Qmax — Q] ( ■ 

[ e=k ) 

At time 0, we have: 

P(0,5o,0) = sup E (Y,e-''^qt,(St, - K)] . 

(."ytfe JO<fc<n-lfc.AQ Q \fc = / 



n-1 



3 



Note that this corresponds to the hmit case of the contract with penahzed constraints when A = 
B = +00. Furthermore, one shows that when the penalties A,B ^ +00 in ([2|), the "penahzed" 
price converge to the "firm" price. This has been confirmed by extensive numerical implementations 
of both methods. In practice when A, B ~ 10 000 both methods become indistinguishable for usual 
values of the volume constraints. 

1.2 Canonical decomposition and normalized contract 

In this section we obtain a decomposition of the payoff of our swing contract (with firm constraints) 
into two parts, one having a closed form expression. It turns out that this simple decomposition 
leads to an impressive increase of the precision of the price computation. It plays the role of a vari- 
ance reducer. Moreover, its straightforward financial interpretation leads to a better understanding 
of the swing contract. 

In fact, we can distinguish a swap part and a normalized swing part: 

P(0,5o) = E(|](7„i„e-'-*'=(5t, -i^)J (4) 

\fc=0 / 



Swap 

/n-1 



+ {qmax - qmin) _ SUp ''^ ^ '"^^Itki^tk " K) 



fefe)e^|^™]'"''''"''"(o,o) Vfc=o 



Normalized Contract 



n-1 



where 

^Q™„,Q™a.(^^g) = {{qt,)k<e<n-uqt, : {^,:Ft„F) ^ [0,1],^%, E [{Qmin - Q) ,Qmax-Q]} 
and 

min ^qmin max f^qmiri)+ 

^min — ) Wmax — ■ y^j 

qmax qmin qmax qmin 

The price models investigated in the following sections define the spot price as a process centered 
around the forward curve, and so ]E(5'f) = -Fo,t is known for every [O,?"]- Thus, the swap part 
has a closed form given by 

n-1 

Swapo = qmin X] (^o,tfe - K). 

i=0 

The adaptation to contracts with penalized constraints is straightforward and amounts to mod- 
ifying the penalization function in an appropriate way. 

1.3 Dynamic programming equation 

In [6], it is shown that, in the penalized problem, optimal consumption is the solution of a dynamic 
programming equation. 

Proposition 1.1. Assume that for some positive constants p and C, the following inequality holds 
for any x > 0, and Qg [n qmin, n qmax] 



\PTix,Q)\ <C{l + xP). 
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Then, there exists an optimal Markovian consumption q*{tk, St^,Qt,,) given by the maximum argu- 
ment in the following dynamic programming equation: 

P{tk,St„Qt,) = max {q{St, - K) + e-''^'>'+^-'MPitk+i, St,^,,Qt, + q)\St,)} , , , 
P{T,St,Qt)=Pt{St,Qt). 

Usually, the function Pj,{x, Q) is given by ([2]). Then, the case with firm constraints corresponds 
to the limit case where Pj,{x,Q) = (-oo)1{^^[q^.^^,q^^^]}. 

When considering a contract with firm constraints, a more operating form (see [4]) can be the 
following 

P{tk, St„Qt,) = ma4q{St,-K) + E(P(ifc+i, St,^^ , Qt^ +q)\St,), (7) 

q £ [^min ) <?max ], Qtu+q^ [(Qmin - {n- k) (Qmax - (n - k)qmin)+]} ■ 



1.4 Bang Bang consumption 

1.4.1 Case with penalties on purchased volumes 

[6] showed the following theoretical result. 

Theorem 1.2. Consider the Problem {1\ and Pj,{x,Q) = —xP{Q), P being a continuously differ- 
entiable function. If the following condition holds 

F{e--'>'{St,- K)+E{e-^^STP'{Q*T)\St„Ql) = 0) = 0, 

the optimal consumption at time tk is necessarily of bang-bang type given by 

q*{tk, St^.Ql^) = Q'maxl{e-'-*fc(Stj^-E-)+E(e~'-l'5TP'(Q?.)|Stj^,Q?J>0} 

+gmml{e-'-tfe {St^^ -K)+E{e-^T g^p, (^Q*^)^St^ .Qt\ )<0} " 

The above assumption seems difficult to check since it involves the unknown optimal consump- 
tion. However, this would be the case provided one shows that the random variable e~''*'= (Sj^. — 
K) + E(e~^^S'TP'(Q5^)|5'tj., <5j^) is absolutely continuous as noticed in [6]. 



1.4.2 Case with firm constraints 

In the companion paper [5] , we establish some properties of the value function of the swing options 
viewed as a function of the global volume constraints {Qmin, Qmax)- Thanks to (j4]) one may assume 
without loss of generality that the contract is normalized, i.e. qmin = and qmax = 1- We consider 
the following value function: 

P{Qmin, Qmax) = SUp E [ |^ e-^''^qt, {St, - K)] 

('?tfe)0<fc<n-lt-4Q_() \fc=0 / 

defined on the unit (upper) simplex {(n, G M^, < u < f < n}. 

Proposition 1.3. The premium function {Qmin, Qmax) ^ P{Qmin,Qmax) is a concave, piecewise 
affine function of the global purchased volume constraints, affine on elementary triangles (m, M) + 
{{u,v),0 < u < V < 1}, {m,M) £ N^, m < M < n and {m,M) + {{u,v),0 < v < u < 1}, 
(m, M) S N^, m<M— l< 71 — I which tile of the unit (upper) simplex. 
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Theorem 1.4. For integral valued global constraints, i.e. {Qmin, Qmax) ^ N^; there always exists a 
bang-bang optimal strategy i.e. the a priori [0, l]-valued optimal purchased quantities ql^ are in fact 
always equal to or 1. 

Remark 1.2. This result can be extended in some way to any couple of global constraints when all 
the payoffs are nonnegative (see [3]). Furthermore, it has nothing to do with the Markov dynamics 
of the underlying asset and holds in a quite general abstract setting. 

An example of the premium function {Qmin, Qmax) ^ P{Qmin, Qmax) is depicted on Figured! 




Figure 1: Value function P{Qmin,Qmax) versus the global constraints 

Now we turn to the problem of the numerical evaluation of such contracts. As announced, we focus 
on an optimal quantization algorithm. 

2 Optimal quantization 

Optimal Quantization [211 [H [2l [3] is a method coming from Signal Processing devised to approxi- 
mate a continuous signal by a discrete one in an optimal way. Originally developed in the 1950's, 
it was introduced as a quadrature formula for numerical integration in the late 1990's, and for 
conditional expectation approximations in the early 2000's, in order to price multi-asset American 
style options. 
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Let X be an M"^- valued random vector defined on a probability space ($7,^, P). Quantization 
consists in studying the best approximation of X by random vectors taking at most fixed values 

Definition 2.1. Let x = G (M'')^. A partition (Ci(x))i=i,...,Ar of is a Voronoi 

tessellation of the A^-quantizer x (or codebook; the term grid being used for {x^, . . . ,x^}) if, for 
every i G {1, . . . , N}, Ci{x) is a Borel set satisfying 

C,{x) C {C G M'^, le - x'i < min 1^ - x^\} 

where | . | denotes the canonical Euclidean norm on W^. 

The nearest neighbour projection on x induced by a Voronoi partition is defined by 

Proj^ : 2/ G M"^ 1-^ X* if yG Ci{x). 
Then, we define an x- quantization of X by 

X^ = Proj,(X). 

The pointwise error induced when replacing X by X^ is given by \X — X^\ = d{X, {x^, . . . , x^}) = 
mini<j<Ar \X — When X has an absolutely continuous distribution, any two x-quantizations 
are P-a.s. equal. 

The quadratic mean quantization error induced by the the A^-tuple x G ffi'^ is defined as the 
quadratic norm of the pointwise error i.e. \\X — X^\\^. 

We briefly recall some classical facts about theoretical and numerical aspects of Optimal Quan- 
tization. For details we refer e.g. to |15l I22j . 



3 1 



- 



• • • • 
• • • • • 



. • • • . • . . ' ' • 

• • • •• • 

• •• ••••• •• 

• • • • • • • . . • . • • 



• • • 



• • • 



Figure 2: Optimal quadratic quantization of the normal distribution AA(0, 12), N = 200. 

Tiieorem 2.1. Let X€ L^(M'^,P). The quadratic quantization error function 
X = (x\...,x^) 1 — >E( min |X-x*|^) = ||X-X^|L 

reaches a minimum at some quantizer x* . Furthermore, if the distribution ¥x has an infinite 
support then x*'^^^ = (x*'^, . . . , x*'^) has pairwise distinct components and N i—>- min^^^^d-jN \\X — 
X^IIj is decreasing to as N ] +00. 
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Figure [2] shows a quadratic optimal quantization grid for a bivariate normal distribution A^(0, 12)- 
The convergence rate to of optimal quantization error is ruled by the so-called Zador Theorem. 

Theorem 2.2. fl^ Let X £ L^+\F),6 > 0, with FxidS,) = ip{^)Xd{dC) + fidS,), ± Xd (\d 
Lebesgue measure on M.'^). Then 



lim (iV^ min \\X - X'^ll^) = J2 d ( [ ^f^+^dXd] 



The true value of J24 is unknown as soon as d > 3. One only knows that J24 = + o{d). 
Zador's Theorem implies that \\X - ' \\^ = 0{N~d) as N ^ +00. 

Proposition 2.3. i21\ \22^ Any L'^-optimal quantizer x £ satisfy the following stationarity 
property 

E{X\X^) =X^. 
In particular, for any stationary quantizer E{X) = E(X^). 

The random vector X^ takes its value in a finite space {x^, . . . ,x^}, so for every continuous 
functional / : M*^ ^ R with f{X) G L^{F) , we have 

AT 



E(/(X-)) = 5^/(x*)P(XgC,(x)) 



which is the quantization based quadrature formula to approximate E(/(X)) |2H I22j. As X^ is 
close to X , it is natural to estimate E(/(X)) by E(/(X^)) when / is continuous. Furthermore, 
when / is smooth enough, on can upper bound the resulting error using ||X — X^Hj, or even 
||X — X^llj (when the quantizer x is stationary). 

The same idea can be used to approximate the conditional expectation K{f(X)\Y) by K(f{X)\Y), 
but one also needs the transition probabilities: 

nX G C,(x)|y G Ci{y)). 

The application of this technique to the quantization of spot price processes is discussed in 
details in the Annex, page[ 



3 Pricing swing contracts with optimal quantization 

3.1 Description of the algorithm (general setting) 

In this section we assume that {St^,)Q<k<n is a Markov process. For sake of simplicity, we consider 
that there is no interest rate. We also consider a normalized contract, as defined in Section [1.2i 

In the penalized problem, the price of the swing option is given by the following dynamic pro- 
gramming equation (see Equation \^ : 

Pit,,St„Qt,) = max [q{St,-K)+E{P{tk+i,St,^„Qt,+q)\St,)] 
ge{o,i} 

P{T,ST,QT)=Pr{ST,QT) 

where tk = kA, k = 0, . . . ,n, A = — . 
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The bang-bang feature of the optimal consumption (see Section ll.4p allows us to limit the 
possible values of q in the dynamic programming equation to q£ {0, 1}. At time tk, possible values 
of the cumulative consumption are 



Q 



,0<l<k. 



At every time we consider a(n optimized) A^^-quantization 5"^,, = S^^ , k = 0, 
on an optimized quantization A'^-tuple (or grid) x^^'""^ := ^s^, . . . , s^*^ of the spot St,^. 

The modeling of the future price by multi-factor Gaussian processes with memory (see Section 4 
for a toy example) implies that (5't)tg[o,r] is itself a Gaussian process. Then the quantization of 5tj. 
can be obtained by a simple dilatation-contraction (by a factor StD(S'tj,)) from optimal quantization 
grids of the (possibly multivariate) normal distribution, to be downloaded on the website |23| 

www. quantize .maths-fi . com 

Then we compute the price at each time tk, for all points on the corresponding grid, and for all 
the possible cumulative consumptions: 



(8) 

, n based 



f P{tk,sl,Qt,)= max [q{si-K)+E{P{tk+i,St,^,,Qt,+q)\St,=si) 
cje{o,i} 

1, 



(9) 



1, 



1 = i,... ,Nk, 
{ P{T,s'j.,Qt) = PA4,Qt), 

When considering a contract with firm constraints, we need to compute the price at each time 
tk, for all the points of the quantization grid of the spot price, and for all the admissible cumulative 
consumptions (See Figure [3j) 



Q 



{n - k + 1)) , ,£ = 0, min(/c, 



(Qr 



{n-k + 1)),. 



(10) 




Admissible 
Volumes 



40 60 
Days 




100 



Figure 3: Volume Constraints 



Using the bang-bang feature (See section [L^ and the dynamic programming principle (Equation [7]), 
this price is given by 



Pitk,sl,Qt,)=ms.x[q{sl-K)+E{P{tk+i,St,^„Qt,+q)\St, = si), (11) 
q£ {0, 1}, Qtk+Q^ [(Qmin - (n - A;))_^ , Qmaxlj • 
Since S't^, takes its values in a finite space, we can rewrite the conditional expectation as: 

Nk+i 

E(P(ifc+i,4^,,Q)|5i, = 4) = J2 



where 

<'' = n4+. = 4+ii4 = 4) 

is the quantized transition probability between times tk and tfc+i- The whole set of quantization 
grids equipped with the transition matrices make up the so-caUed "quantization tree". The tran- 
sition weights (vr^"') matrices are the second quantity needed to process the quantized dynamic 
programming principle ([Ul or [TT|) . A specific fast parallel quantization procedure has been devel- 
oped in our multi-factor Gaussian framework to speed up (and parallelize) the computation of these 
weights (see Annex) . In a more general framework, one can follow the usual Monte Carlo approach 
described in [1] to compute the grids and the transition of the global quantization tree. 

When (St^) has no longer a Markov dynamics but appears as a function of a Markov chain 
{Xk): 5tj, = f{Xk), one can proceed as above, except that one has to quantize (X^). Then the 
dimension of the problem (in term of quantization) is that of the structure process (Xk). 

Of course one can always implement the above procedure formally regardless of the intrinsic 
dynamics of (5"^^). This yields to a drastic dimension reduction (from that of X downto that of 
(St^)). Doing so, we cannot apply the convergence theorem (see Section [3?3|) which says that in a 
Markovian framework the premium resulting from ([9]) or (jlip will converge toward the true one as 
the size of the quantization grid goes to infinity. 

This introduces a methodological residual error that can be compared to that appearing in [5] 
algorithm for American option pricing. However, one checks on simulations that this residual error 
turns out to be often negligible (see Section 14. 3p . 

3.2 Complexity 

The first part of the algorithm consists in designing the quantization tree and the corresponding 
weights. The complexity of this step is directly connected to the size of the quantization grids 
chosen for the transitions computation in 1-dimension, or to the number of Monte Carlo simulations 
otherwise. However those probabilities have to be calculated once for a given price model, and 
then several contracts can be priced on the same quantization tree. So we will mainly focus on the 
complexity of the pricing part. 

We consider a penalized normalized contract, i.e Qmin = and Qmax = 1- The implementation 
of the dynamic programming principle requires three interlinked loops. For each time step k (going 
backward from n to 0), one needs to compute for all the points s\, i = 1,. . . ,Nk of the grid and 
for every possible cumulative consumption Q^^ (0 < ^ < /c) (see dH)) the functional 

max [q{sl - K) + E(P(tfc+i, 5^,^, , Qt, + q)\St, = si)] 

(?G{0,1} 

which means computing twice a sum of N^+i terms. 
Hence, the complexity is proportional to 

n-l 

In the case where all layers in the quantization tree have the same size, i.e N = Nk,'\/k = 1, . . . , n, 
the complexity is proportional to " ^ . This is not an optimal design but only one grid needs to 
be stored. It is possible to reduce the algorithm complexity by optimizing the grid sizes Ni ^(with 
the constraint iV^. = nN), but it costs more memory space. 



To minimize the complexity, set Nk ~ (k+^"iog(n) ' = 0, . . . , n— 1, which leads to a global complexity proportional 
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In the case of firm constraints, the dynamic programming principle (jlip has to be computed 
for every admissible cumulative consumption, i.e for every Ql^ (0 < £ < min(A; — l,Qmax) — 
(Qmin — {n — k + !))_,_, see (fTUj) ). The complexity is proportional to 

n-l 

J^(min(fc, g^ax) - (Qmin - {n - k + 1))+ + l)NkNk+l. 
k=0 

The complexity in the case of firm constraints is lower than the one for a penalized problem, 
and depends on the global constraints (Qmim Qmax)- But the implementation is easier in the case 
of a penalized problem, because one does not need to check if the cumulative consumption volume 
is admissible. Both approaches have been numerically tested and results are indistinguishable for 
large enough penalties. For the implementation readiness, the approach with penalties has been 
adopted. 

In order to reduce the complexity of the algorithm, one usually prunes the quantization tree. 
In most examples, at each layer k, many terms of the transition matrix are equal to 

or negligible. So while the transition probabilities are estimated, all the transitions that are not 
visited are deleted. This step is important because it allows to reduce significantly the algorithm 
complexity. 

In practice we can even neglect transitions whose probability is very low, say less than 10~^. 



3.3 Convergence 

In [3] is proved an error bound for the pricing of swing options by optimal quantization. 

Let Pq{Q) denote the price of the swing contract at time 0. n is the number of time step, 
and Q = {Qmin,Qmax) IS the global constraint. We consider a contract with normalized local 
constraints, i.e qmin = and Qmax = 1- The "quantized" price Pq{Q) is the approximation of the 
price obtained using optimal quantization. 

Proposition 3.1. Assume there is a real exponent p G [l,+oo) such that the (d- dimensional) 
Markov structure process {Xk)o<k<n-i satisfies 

max G LP+^(P),r? > 0. 

0<A;<n— 1 

At each time k G {0, . . . ,n — 1}, we implement a (quadratic) optimal quantization grid of size 
NofXk. Then 

71 

II sup \p-{Q)-P^{Q)\\\p<C— 

QeT+{n) Nd 

where T'^{n) := {(ii,i)),0 <u<v<n} is the set of admissible global constraints (at time 0). 

In fact this error bound turns out to be conservative and several numerical experiments, as 
those presented in Section HI suggest that in fact the true rate (for a fixed number n of purchase 
instants) behaves like 0{N~d). 



4 Numerical experiments 

In this section the same grid size has been used at each time step, i.e. we always have = N,k = 
0, . . . , n. The results have been obtained by implementing the penalized problem and using the 
canonical decomposition (see Section [L2]) . 
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4.1 The one factor model 



Swing options are often priced using the least squares regression method "a la Longstaff-Schwartz" 
|19j . This section aims to compare our numerical results to those obtained with LongstafF-Schwartz 
method. We consider a one factor model, which corresponds to a one dimensional Markov structure 
process. 



4.1.1 Quantization tree for a one dimensional structure process 

We consider the following diffusion model for the forward contracts iFf x)o<t<T- 



Ft,T 

where Vl^ is a standard Brownian motion. It yields: 



St = Fo,t exp (^a e""(*-^)(iTy, - 



where 



2 

2 /I „-2at\ 



Denote Xt = /(, e-"(*-")dW,. The structure process {Xk)k>Q can be quantized using the fast 

parallel quantization method described in the Annex (page [22]) . Let x^^"^ denote an (optimal) 
quantization grid of of size A^. We have to compute for every k E {0, . . . ,n — 1}, and every 
(i, j) S {1, . . . , A^}^, the following (quantized transition) probabilities: 

= F{m G C,(xf )); ak+m + G C,{xi1\)) (12) 

where (r?i,?72) ~ M{0,l2), and Ok and (3^ are scalar coefficients that can be explicited. 

This can be done by using quantization again and importance sampling as presented in the 
Annex, page [22] (see (j2T]) ). 



4.1.2 Comparison with the regression method 

We first use the following parameters for the one factor model: 

a = 70%, a = 4, Fo,i, =20, k = 0,...,n. 

The following tables present the results obtained with Longstaff Schwartz and optimal quan- 
tization, for different strike values. 1000 Monte-Carlo sample paths have been used for Longstaff- 
Schwartz method and the confidence interval of the Monte Carlo estimate is given in the table. 
A 100-point grid has been used to quantize the spot price process, and the transitions have been 
computed with a 500-point grid. The local volume constraints qmin and Qmax are set to and 6 
respectively . 

We first consider a case without constraints (Table [1]), which means that the swing option is a 
strip of calls, whose price is easily computed with the Black-Scholes formula. 



Table [2] presents the results obtained with the global constraints Qmin = 1300 and Qmax = 1900. 
Volume constraints are presented on Figure [4] 
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K = h 


= 10 


= 15 


K = 20 


Longstaff-Schwartz 


[32424,33154] 


[21360,22127] 


[11110,11824] 


[3653,4109] 


10 point grid 


32726 


21806 


11311 


3905 


20 point grid 


32751 


21834 


11367 


3943 


50 point grid 


32759 


21843 


11380 


3964 


100 point grid 


32759 


21843 


11380 


3964 


200 point grid 


32761 


21845 


11382 


3967 


Theoretical price 


32760 


21844 


11381 


3966 



Table 1: Comparison for a call strip (no global constraints) 
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Figure 4: Consumption constraints 





K = h 


K = IQ 


K = 15 


K = 20 


Longst aff- S chwar t z 


[29068,;29758] 


[19318,;19993] 


[10265,;10892] 


[2482,;3038] 


10 point grid 


29696 


20216 


10981 


3067 


20 point grid 


29494 


20018 


10841 


2863 


50 point grid 


29372 


19895 


10729 


2718 


100 point grid 


29348 


19872 


10704 


2687 


200 point grid 


29342 


19866 


10698 


2680 



Table 2: Comparison with constraints 



The results seem consistent for both methods, the price given by quantization always belongs 
to the confidence interval of the Longstaff-Schwartz method. One can note that it is true even for 
small grids, which means that quantization gives quickly a good price approximation. Moreover, 
the price given by quantization is very close of the theoretical price in the case of a call strip. 
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4.1.3 Execution time 

In this section are compared the execution times to price swing options using optimal quantization 
and Longstaff-Schwartz method. 

The size of the quantization grid is 100 for the pricing part and 200 for the transitions compu- 
tation. And 1000 Monte Carlo simulations are used. The maturity of the contract is one year. 

The computer that has been used has the following characteristics: 

Processor: Celeron; CPU 2,4 Ghz; 1,5 Go of RAM; Microsoft Windows 2000. 

The execution times given in Table [3] concern the pricing of one contract, which yields the 
building of the quantization tree and the pricing using dynamic programming for quantization. 



Longstaff-Schwartz 


Quantization: 
Quantization tree building -|- Pricing 


Quantization: 
Pricing only 


160 s 


65 s 


5 s 



Table 3: Execution time for the pricing of one contract 



If we consider the pricing of several contracts, there is no need for re computing the quantization 
tree if the underlying price model has not changed. That is why quantization is really faster than 
Longstaff-Schwartz in this case, as one can note from the results presented in Table HI 



Longstaff-Schwartz 


Quantization 


1600 s 


110 s 



Table 4: Execution time for the pricing of 10 contracts 



4.1.4 Sensitivity Analysis 

When contracts such as swing options ought to be signed, negotiations usually concern the volume 
constraints. That is why the valuation technique has to be very sensitive and coherent to constraints 
variation. In this section we will compare the sensibility to global constraints for Longstaff-Schwartz 
method and optimal quantization. 

Figure [5] represents the price of the contract with regards to the global constraints Qmin and 
Qmaxi and Figure E] represents the price versus Qmax for a fixed value of Qmin equal to 1300. 



One can notice that the surface obtained with optimal quantization is very smooth. If Qrnax 
increases, the price increases. However, it is not always true with Longstaff-Schwartz because of 
the randomness of the method, and the limited number of Monte Carlo simulations imposed by 
the dimension of the problem and the number of time steps. 

New Monte Carlo simulations are done for each different contract, i.e each time Qmin or Qmax 
varies. Of course, the same simulations could be used to price all the contracts, but unfortunately 
these simulations could be concentrated in the distribution queues and give a price far from the 
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Optimal Quantization 
Longstaff-Schwartz 




Figure 5: Sensitivity to global constraints 



Optimal Quantization 
Longstaff-Schwartz 




Fi gure 6: Sensibility to Q^i 



real one for all the contracts. As concerns quantization, the grid is build in order to give a good 
representation of the considered random variable. One of the great advantages of optimal quanti- 
zation over Monte-Carlo is that this first algorithm always approximates the whole distribution of 
the payoff meanwhile it can take a while before Monte-Carlo explores some parts of it. 
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4.1.5 Convergence 



In this section we will study the convergence of the quantization method. We focus on the conver- 
gence of the pricing part of the algorithm. 

We consider a one year maturity contract with the volume constraints depicted on Figure HI 
and the daily forward curve depicted on Figure [71 



35 
30 




5 

J ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ r- 

oooooooooooo 
oooooooooooo 

OOOOOt-OOOOOt- 



Figure 7: Daily forward curve 



Let P{N) be the price obtained for a quantization grid of size A^, the error has been computed 
as \P{N) — P(400)|, N < 400. We assume that the error can be written as a functional of the grid 
size N with the following shape: 

A linear regression in a logarithmic scale is done to find the functional that best fits the empirical 
error. The a coefficient obtained is 1.96. 

Figures [8] and [9] show the obtained numerical convergence and the corresponding fitted func- 
tional. 
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Figure 8: Numerical Convergence 
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Figure 9: Numerical Convergence (logarithmic scale) 

The same experiments have been done for other contracts, results are presented in Table [5j 
qmin and qmax are set to and 6. 



Forward Curve 


Strike 


Constraints {Qmin - Qmax) 


Estimated a 


Figure [7] 


20 


1300-1900 


1.96 


Flat (20) 


20 


1300-1900 


2.07 


Flat (20) 


10 


1300-1900 


2.32 


Flat (20) 


20 


1000-2000 


1.95 


Flat (20) 


20 


1600-1800 


2.26 



Table 5: Estimation of the convergence rate 



We can conclude that the convergence rate of the quantization algorithm for pricing swing 
options is close to 0{^). This convergence rate is much better than Monte-Carlo, and leads to 
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think that optimal quantization is an efficient alternative to Longstaff-Schwartz method for this 
problem. 



4.2 Two factor model 

We consider the following diffusion model for the forward contracts {Ft^T)o<t<T- 



dFt,T 

Ft,T 



where and are two Brownian motions with correlation coefficient p. 
Standard computations based on Ito formula yield 

St = Fo,t exp (^ai e^'^'^^-'Uwl + ^2 ^ e~'^^'^*~'Uw'^ - ^A^^ 

where 

A? = ^(1 - e-2"i*) + ^(1 - e-2"2*) + 2p^i^(l - e-("i+"2)*). 
2ai 2a2 ai + 02 

Unlike the one factor model, the spot price process obtained from the two factor model is not a 
Markov process. Hence the dynamic programming equation ^ cannot be used directly. However, 
the structure process of the two factor model (See Annex, page [25]) 



is a Markov process, and St = f{Xt) where / : 1-^ M is a continuous function. So we can rewrite 
the dynamic programming equation as follows: 

P{tk,Xt,,Qt,) = max,g[,_,,_] {q{f{Xt,) - K) +nP{tk+i,Xt,^,,Qt, + q)\Xt,)], 

(13) 

P{T,XT,QT)=Pr{f{XT),QT)- 

Then we need to quantize the valued structure process (Xfc)fc>o. This can be done using 
the Fast Parallel Quantization (See Annex). The transitions are computed using Monte Carlo 
simulations and importance sampling (see p^ ). 

4.2.1 Call strip 

We first consider a case without constraints, and compare the results with the theoretical price 
of the call strip, for several values of the strike K. The maturity of the contract is one month. 
Transitions have been computed with 3000000 of Monte-Carlo simulation. The parameters of the 
two factor model are: 

0-1 = 36%, ai = 0.21, CJ2 = 111%, 02 = 5.4, p = -0.11. (14) 

Table [6] presents the results obtained for a strip of call. Even if the quantized process is taking 
values in M^, prices are close to the theoretical price even for small grids. 
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K = b 


= 10 


= 15 


K = 20 


Theoretical Price 


2700 


1800.21 


924.46 


268.59 


50 point grid 


2695.26 


1795.26 


918.18 


261.06 


100 point grid 


2699.79 


1799.89 


923.64 


267.17 


200 point grid 


2697.67 


1797.83 


921.94 


266.71 


300 point grid 


2702.02 


1800.16 


924.40 


268.52 



Table 6: Call Strip 
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Figure 10: Numerical Convergence 



Forward Curve 


Strike 


Constraints {Qmin - Qmax) 


Estimated a 


Flat (20) 


10 


80-140 


1.26 


Flat (20) 


20 


80-140 


1.00 


Flat (20) 


20 


30-170 


0.67 


Flat (20) 


20 


100-120 


1.19 



Table 7: Estimation of the convergence rate 



4.2.2 Convergence 

We use the same procedure as in section 14.1.51 to find the functional N ^ that best fits the 
empirical error. 

Figure [10] shows an example of the empirical error and table [7] gather the values of a obtained 
for different contracts. The contract maturity has been set to one month. 



The convergence of the quantization algorithm is close to 0{jj). The convergence rate is linked 
to the dimension d of the structure process, and from the results obtained in section 14.1.51 and in 
this section, we can assume that the convergence rate is close to 0( jy2/d )i which is better than the 
error bound theoretically established in Section [3.31 
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4.3 Dimension reduction 

In the case of multi- factorial models, we need to quantize the structure process [X^)]^ instead of 
the spot process 5ij. = /(Xf.) in order to work with a Markov process (See section Br2]) . Only the 
one factor model is Markovian. That is why quantization and Longstaff-Schwartz method have 
been compared just for this model. Longstaff Schwartz method also requires a Markov underlying 
process. 

From an operational point of view, it is interesting to study the results obtained by formally 
quantizing the spot process (St,.), regardless to its dynamics, and using the approximation 

EiX\TtJ^E{X\St,) 

for any random variable X, even if the spot price (St^) is not a Markov process. Similar approx- 
imation has already been proposed by Barraquand-Martineau in [5]. Numerical tests have shown 
that the resulting prices remain very close to those obtained by quantizing the structure process in 
the case of a two factor model. Execution time and convergence rate are significantly faster, and 
the quantization tree can be computed as presented in the Annex page 1221 using Equation (j2ip . 

Even if there is no theoretical evidence on the error, this approach seems useful to get quick 
results. Table [8] presents some results, parameters of the two factor model are those of ()14p and 
volume constraints are represented on Figure HI 



Quantized Process 


K = 5 


K = 10 


K = 15 


K = 20 


Spot Price 


30823.38 


21414.30 


13021.87 


5563.68 


Structure Process 


30705.78 


21518.75 


13123.56 


5722.89 



Table 8: Quantization of the spot price vs Quantization of the bivariate structure process 

The convergence rate obtained in this case is always 0(-^), which is consistent with the general 
rate 0( ^2/d )i because the quantized process is M- valued. 

Therefore, even if the spot process is not Markov, the quantization method can be performed 
all the same way as if it were, with small damage in practice. This is of course not a general 
consideration but rather an observation over the considered problem. The dramatic increase in the 
computation effort that can be gained from this observation can justify in this case a lack of formal 
rigor. 

Conclusion 

In this article, we have introduced an optimal quantization method for the valuation of swing op- 
tions. These options are of great interest in numerous modeling issues of the energy markets and 
their accurate pricing is a real challenge. 

Our method has been compared to the famous Longstaff-Schwartz algorithm and seems to perform 
much better on various examples. In fact, the optimal quantization method shares the good prop- 
erties of the so called tree method but is not limited by the dimension of the underlying. Moreover, 
specific theoretical results provide a priori estimates on the error bound of the method. 
Thus, optimal quantization methods suit very well to the valuation of complex derivatives and 
further studies should be done in order to extend the present results to other structured products 
arising in the energy sector. 
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Annex: Fast Parallel quantization (FPQ) 

In this annex we propose an efficient method to quantize a wide family of spot price dynamics 
(5fj.). To be precise we will assume that a time discretization A being fixed, 

Skii = fk{Xk),k>{) (15) 

where {Xk)k>o is a M"*-valued Gaussian auto-regressive process and {fk)o<k<n a family of contin- 
uous functions. The fast quantization method applies to the Gaussian process (^A;)o<A;<n- We will 
apply it to a scalar two factor model in full details. As a conclusion to this section we will sketch 
the approach to a multi-factor model. 

Quantization of the Gaussian structure process 

We consider a centered Gaussian first order auto-regressive process in R"^: 

Xk+i = AXk+Tek+i (16) 

where A G M{m x m,M), T G A4{m x m,R) lower triangular, and (sfc) i.i.d. with N{Q,Im) distri- 
bution. 

Denote by D{Z) = [K{ZiZj)]-^^ - the covariance matrix of Z. We have Vfc G N: 

D{Xk+i) = AD{Xk)A* +TT*. 
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Denote the lower triangular matrix such that D{Xk) = S^E^. 

We consider for every /c = 0, . . . , n — 1, an optimal (quadratic) quantizer x^^^^ of size N^, for 
the AA(0, Im) distribution. The quantization grid of the random variable is taken as a dilatation 
of 4^''\i.e 

To calculate the conditional expectations in the dynamic programming equation, we need to 
get the following transition probabilities: 

= nXk+i G C^{xk+i)\Xk G Ci{xk)) 
where Ci{x) denotes the i-th Voronoi cell of the generic quantizer xG (M'^)^. Then 

¥{Xk G Ci{xk)) = nZ G C,(xf *^)), 

with Z ~ J\f{0,lm)- This probability is provided as a companion parameter with the normal 
distribution grid files (available on [23]). 

To get the transition probability vr^-' we need to compute 

pt^' = F{Xk+i G Cj{xk+i),Xk G Q{xk)). 

Proposition 4.1. Let X be a discrete time process described as above. Let U , V be two gaussian 
random variables M{0, /m)- Then we have for every k G {0, . . . ,n— 1}, every i G {1, . . . , N^}, every 
j G {!,..., A^fc+i}, 

P(Xfc+i G Cj{xk+i),Xk G Ci{xk))=F{U G Ci{x\!'''^),Ak+iU + Bk+iV G )) (17) 

where Aj. and are q x q matrices whose coefficients depend on k, and on the matrices A, T. If 
k = 0, 

P(Xi G C,{xi)) = F{V G C,{xf (18) 

Proof. We have: 

P(Xfc+i G C,ixk+i);Xk G aixk)) = W{AXk+Tek+i G Cj{xk+i);Xk G Q(xfc)). 

We consider the couple (X^ , Te/c_|_i ) . Te^+i is independent of X^. Let ry = {r]i,r]2) a couple of 
independent Gaussian random vectors: rji ~ N{0,lm),i = 1,2. Then (Xfc,Tefc_|_i) ~ {'Ekr]i,TT]2) 
and 

P(Xfc+i G Cj{xk+iy,Xk G Ci(a;fc)) = P(ASfcr/i + Tr]2 G Cj(xfc+i); Sfc7?i G ^(xfc)) 

= P(S-^, (AS,r?i + r^2) G C,(xi^\+^));77i G C.(xf '=))). 

Setting 
we get 

¥{Xk+i G C,(rEfc+i);Xfc G C,(xfc)) = n^k+m + Bk+m G C,-(xi'^'i+^^); ryi G Ci(xf'=^)). 
If fc = and So = 0, the quantity 

P(XfcGC,-(xfc))=P(7?iGC,-(xf^^)) 
is given as a companion parameter with the quantization grids of the normal distribution. □ 

Remark 4.1. Equation (|17|) emphasizes the fact that the transitions can be computed in parallel. 

Remark 4.2. To simplify the structure of the quantization tree we propose to consider the same 
normalized grid of size A^^ = iV at each step k but other choices are possible like those recommended 
in[l]. 
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Numerical methods 



Hereafter we will focus on the numerical computation of these transitions. 

• The standard Monte Carlo approach The simplest way is to use a Monte Carlo method. 
One just needs to simulate couples of independent gaussian random variables (r/i, 1^2). This approach 
can be used whatever the dimension m of the random variables iji and 772 is. It can clearly be 
parallelized as any MC simulation but fail to estimate the transition form states which are not 
often visited. 

• Fast Parallel Quantization Method In order to improve the accuracy, especially for the 
points of the grids x^j!^^^^ which are rarely reached by the paths starting from the cell of 
2;(A^fe).*^ it is possible to perform importance sampling. The idea is to use Cameron-Martin formula 
to re-center the simulation: for every k G {0, . . . , n — 1} and every i£ {!,..., Ni^}, 

Then these expectations can be computed by Monte Carlo simulations, the transitions between the 
different times steps can be computed in parallel. 

• Quantized Parallel Quantization Method If m = 1, the transitions can be computed 
using again optimal quantization, because in low dimension (say < 4), quantization converges 
faster than Monte Carlo method. In this case, we have to compute a two dimensional expectation. 

We estimate for every k£ {0, . . . , n — 1}, every i G {1, . . . , N^} and every j G {1, . . . , Nf^^i} the 
following probabilities: 

= p(^^ e ^.(xf^)); ak+m + Pk+iV2 G C,{xifr^)) (20) 
where (771, r/2) ~ AA(0, 12), and Ok and f3k are scalar coefficients satisfying af, + (3l = 1: 



k + l 

T 



S 



k + l 



To alleviate notations, we temporarily set x = x^^*"^ and y = x^j!^-^^\ 



X 2,X^2 



[^(x* + 2;* i(x* + x*"*"-^)] = Ci{x) (the same shortcut is implicitly 



We define 
defined for y). 

In order to reduce the problem dimension, it is possible to write the probability double 
integral, and to integrate first with respect to the second variable by using Fubini theorem: 



G Ci{x); ak+im + (3k+im G Cj{y)) 



i-- i+- 



Ell .^1 ^l™i)_AA(i^- '^^±121) 



where N{x) is the distribution function of the normal distribution. 
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Importance sampling can again be used to improve the results precision. Eventually we have 
to compute the following one dimensional expectation: 



'(r/i G Ci{x); ak+im + (3k+im G Cj{y)) 
-(^')^ „ 



ctk+iivi + 



(21) 



For this one-dimensional expectation computation, quantization can be used again since it converges 
faster than Monte Carlo method. 

Example: Two factor model 

We consider the following diffusion model for the forward contracts {Ft^T)o<t<T- 

dFtT 



tT 



where and W"^ are two Brownian motions with correlation coefficient p. 
Standard computations based on Ito formula yield 



St = Fo,t exp (^ai e-'^^'^'-^Uw} + a2 



-"2(t-^)^w2 



-A 



where 



2a2t\ 



+ 2p- 



(J1CJ2 



.(1-e- 



2ai 2a2 ' ' ' ai + 02 

We have St = Fo^texp (cxiX/ + a2Xf — ^A^), where Xt is the following structure process: 



Xt 



(22) 



Proposition 4.2. Let Z = (Zt) be an Ornstein-Uhlenbeck process. Zt = Zo + /(Je-"(*-")(iB, where 
B is a standard Brownian motion, and Zq is Gaussian and independent of B. Z can be written at 
discrete times kA as a first order auto-regressive process: 



Zk+i = e-'^^Yk + Vl-e- 



2aA 




(23) 



where {sk) is i.i.d and ei ^ M{0, 1). 

Xt is made up with two Ornstein-Uhlenbeck processes. Using Proposition 14. 2^ it yields: 
Proposition 4.3. 



Xk+i — AXk + Tsk+i 



with Ek ^ M{0, 12) i.i.d. and 



-aiA 



A 



-02A 
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1 (-1 _ p-(ai+a2)A\ 

_ 01+02^ I 



4ai 

Hence it is possible to use the fast parallel quantization method described in section [43l 
General multi-factor Gaussian model 

More generally, we consider a family of price dynamics that can be written as follows: 

dFtT 



^P,(T-t)e-"»(^-*)dW/ (24) 



Ft.T 

where Pi{x) is a polynomial function of degree dj, for every i = 1,. . . ,m, and is a Brownian 
motion, with d < W\ >t= Pijdt. 

The two factor model (Section I4.3P corresponds to m = 2, Pj = ctj, i = 1,2. 

In order to price a swing option with such a model, we first need to quantize it. Equation (|24p 
yields: 

F,j, = F^^Te^T=.!!,P^{T-s)e-.i^-^Uwi-\4>{t,T) (25) 

where 

0(t,r)=V / PfiT - s)e-'^''^^'^-'^ds + y2p,j / Pi(r-s)P,(T-s)e-("»+°^)('^-^)ds. 
i=i -^0 i^j -^0 

Practically we focus on the spot price Ft^t or the day-ahead contract Ft,t+i- Unfortunately these 
processes are not Markovian in a general setting, except when m = 1 and di = (Ornstein- 
Uhlenbeck process). 

We consider a discretization time step A > 0, and we set, for all i G {1, . . . , m} and for all 
I G {0,...,di} 

f-fcA 



X'/= / Pi((fc + OA-s)e-"'«^'+')-^)dWi. 
Jo 



10 

Proposition 4.4. Xk = [Xfc']i<i<m,o<i<di is a R'^^^""^'^"^'^'^ -valued gaussian AR(1). 

Lemma 4.5. Let P G d^P = d andOeR*. Then {P{Z + /6l))o<Kd a basis ofR'^. 

Proof. Using a dimension argument, only the linear independence of the family has to be checked. 
And we have 

Y,x,Piz+k9)=o ^ x:^^-E^^^'^(^)=o 

k=0 k=0 j=0 

d „j / d \ 
u-' 



Since (p(j)(Z))p< .^^ is a basis of [Z] , it yields 

d 

ViG{0,...,4, J2Xkk^ = 

k=0 

so that Afc = 0, < k < d since det[/c-^]o<A:,j<d / (Vandermonde determinant). □ 
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Proof, (of Proposition 14. 4p We can extend the definition of X^,'' to ^ G N. It is easy to check that 

r{k+l)A 

Jo 



X 



i,l 
k+1 



Pi iik + l + I) A - s) e-°«{{^+i+OA-s)^p^^i 



-5^i,«+l , i,l 



where e^^^^ = if^^^^ Pi ((fc + 1 + /)A - s) e-^^i.ik+l+l)A~s) ^^i^ 
If I = di, 



According to Lemma l4.5t 



\ri,di -irijdi+l . i,di 



P, {Z + {d, + 1)A) = X'''Pi{Z + /A). 

1=0 



Hence 



X 



i,di+l 



kA 



Pi {{k + d^ + 1)A - s)e- 



-ai{{k+di+l) As) 



VA^-'/ Pi{{k + l)A 

1=0 

di 

^i,l^-ai(di+l-l)Aj^i,l 

1=0 
di 

E^'^k'- 



s e 



,{{k+l) As) ^^r^^-a,{d,+l-l) A 



1=0 



FinaUy we have 
where 



■^k+l — ^^^k + ^k+1 

/ 1 ••• \ 



A' 









1 
. .. ~xiA J 



and 



^k+i — + ek,Xo — 

where Sk & a {W^ — Wl^, kA < u < {k + 1)A, z = 1, . . . , m) is independent of J-'^- The process 
{Xk)k is thus a gaussian AR(1). □ 

Xk is the structure process for the spot price SkA = FkA,kA- Its dimension is i(c^j + !)• For 
the two factor model, the structure process is M^-valued, because m = 2, and = 0, i = 1, 2. This 
is coherent with (1221). 
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